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Fluctuation theorems have a very special place in the study of non equilibrium dynamics of physical systems. 
The form in which it is used most extensively is the Gallavoti-Cohen Fluctuation Theorem) 1] which is in terms 
of the distribution of the work p(W)/p{—W) = exp(a?VF). We derive the general form of the fluctuation 
theorems for an arbitrary Gaussian Markov process and find conditions when the parameter a becomes a uni¬ 
versal parameter 1/kT. As an application we consider fluctuation theorems for classical cyclotron motion of an 
electron in a parabolic potential. The motion of the electron is described by four coupled Langevin equations 
and thus is non-trivial. The generalized theorems are equally valid for non-equilibrium steady states. 

PACS numbers: 05.40.-a, 05.70.Ln, 05.60.Cd 


I. INTRODUCTION 

Recent years have seen an upsurge of interest in general¬ 
ized fluctuation dissipation like theorems for driven, nonequi¬ 
librium, dissipative systems[2-8]. Much of this renewed in¬ 
terest has concentrated on a new time dependent dynamical 
variable W (t) which characterizes the work done on the sys¬ 
tem by an external, time dependent force. The natural, con¬ 
jugate (to the applied force) variable is the velocity or the 
current. We re-examine the issue of the fluctuation theorem 
(FT) when the velocity and force fields are multidimensional. 
Typically the stochastic dynamics of many systems can be ad¬ 
equately described in terms of a Gaussian Markov process as 
for example is the case of the Brownian motion. After focus¬ 
ing on some general results we present specific, illustrative 
expressions for the much studied problem of dissipative cy¬ 
clotron motion, which is necessarily two-dimensional in the 
plane normal to the applied magnetic field. Our aim is to criti¬ 
cally assess the validity of the oft quoted fluctuation theorems 
when the two dimensional cyclotron motion is both symmetric 
and asymmetric. 

This FT has been subject to much recent scrutiny, especially 
in the context of more than one-dimension[9-ll]. The ques¬ 
tion asked is: Are fluctuations symmetric? For instance, in the 
case of non-equilibrium fluctuations in two-dimensional sys¬ 
tems, does the GCFT relation apply for forces applied in x and 
y-directions? With regard to this issue, new fluctuation theo¬ 
rems going beyond GCFT, have been proposed that deal with 
so-called ‘isometric’ fluctuations. Recent experiments on a ta¬ 
pered rod undergoing Brownian motion in a sea of spherical 
beads have seen evidence for such isometric fluctuations, due 
apparently to the inherent anisotropy of the system) 1 1]. 

The organization of this paper is as follows-In Sec II we 
introduce the FT for a general N-degree of freedom, classical 
dissipative system, in terms of the work done on the system 
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by an applied time dependent force. We prove that the distri¬ 
bution of work follows a Gaussian as long as the underlying 
stochastic process for the systems dynamics is Gaussian and 
Markovian. In Sec III we derive conditions for the determi¬ 
nation of the parameters that characterize the Gaussian dis¬ 
tribution of the work. We also seek conditions under which 
the parameter a in the FT p[VF]/p[—FF] = exp[akF] has a 
universal form 1/KT. In Sec IV we illustrate our results of 
Secs II and III for the case of the dissipative cyclotron motion 
in a confined parabolic potential. We present explicit forms 
for {W) and fluctuation in W in terms of the temperature T 
and the cyclotron frequency ujc- The results of Sec III, though 
illustrated in Sec IV for the symmetric two dimensional cy¬ 
clotron motion, have much more general validity as discussed 
in the concluding Sec V. 


II. THE FLUCTUATION THEOREM FOR A SYSTEM 
WITH N DEGREES OF FREEDOM 

In this section we give a derivation of the FT for a system 
with N degrees of freedom. We will assume that the system’s 
dynamics is described by a Gaussian Markov process. Let us 

denote the coordinates by xi,X 2 ,X 3 , . xn and kinematic 

momenta, scaled by mass, by vi,V 2 ,V 3 , ....vjy. The Langevin 
equations describing system’s dynamics are assumed to have 
the form 

= Vi, ( 1 ) 

Vi — ^ ^ ^ ^ ^tjVj “f “f Vi • (^) 

3 3 

Here Ay and By contain both dissipative as well as the coher¬ 
ent aspects of the dynamics. The are the delta correlated 
Gaussian random process with zero mean. The fi {t) is the ex¬ 
ternal force that acts on the variable Vi. As an example a force 
on the xi degree of freedom would appear as a driving term 
in the equation for the momentum vi. The work done by the 
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external force on the system is 


B. Determination of a using dynamical equations (1) and (2) 


W = [ dTfj{T)xj{T)mj 

= y] / dTfj{T)vj{T)mj . 

do 


We write Eqs.(l) and (2) using matrices. We define column 
marices ip, N, F with 2N components arranged as ipi = Xi,\ 
1p^+N = Vi,-, N,+n = Fi+n = fi, otherwise Ni and F, 
are zero. Then we write Eqs.(l) and (2) as 

iP = MiP + N + F,. (11) 


A. Distribution of work 


We first show that the distribution of W is Gaussian. To do 
this we construct the characteristic function of W 


C{h) = 

= PpdTfi(T)xi{T)\^ 

— gihiW) fp dTfi{T){vi{T)-(vi{T)))'^ ^ 


where 


{w) = y m. 


drh{j){v^{T)). 


(5) 


It should be borne in mind that {vipr)) depends linearly on 
the external forces. The average appearing in Eq.(4) is just the 
functional for Gaussian process which is well known. Using 
this Eq.(4) becomes 

C{h) = , (6) 

where 

^ dT2fi{Ti)fj{T2)mimj 

Jo Jo 

X ((t^i(n)Wi(T2)) - ('Ui(Tl))(Uj(r2))) . (7) 


HencQp{W), whichis the Fourier transform of C{h), is Gaus¬ 
sian: 


p{W) 


1 

\f^a 


exp{ 


2a2 ^ 


Therefore, 


where 


p{-W) 

p{W) 


exp(—alU), 


a = 2{W)la'^ . 


( 8 ) 


(9) 

( 10 ) 


The Eq.(9) is the generalized Eluctuation Theorem valid for 
any Gaussian Markov process. The parameter a can in prin¬ 
ciple depend on f{t). In many cases discussed in literature 
the parameter a is a universal constant 1 jKT where K is the 
Boltzmann constant and T the temperature. We then recover 
the ET of Gallavoti and Cohen. We will see next that in order 
to obtain the universal constant, the system’s dynamics has to 
satisfy certain specific conditions. 


The matrix M is easily constructed in terms of A and B ma¬ 
trices. Then 

ip{t) = V{t)ip(0)+ f V{t — T)F{T)dT 
Jo 

+ [ V{t-T)N{T)dT, V{t) = e^\ (12) 

and hence 

(^(t)) = f V{t- r)F{r)dr . (13) 

^0 

Using Eq.(l 1) in Eq.(3) we get the mean value of the work 

W = yy / dri [ dT2fi{Ti) 

i J Jo Jo 

X Vi+NJ+Niri - T2)fj{T2)mi . (14) 


Using Eq.(3) we obtain the fluctuation in W as 


a 


2 



dTidT2m,mjfi{Ti)fj{T2) 


X {Vi{Ti)Vj{T2)) , 


(15) 


where the average in Eq.(15) is to be calculated using the 
Langevin equations (11) with U = 0. It can be shown that 


2N 

{v^{Tl)Vj{r2)) = “ T2)\i+nA^I‘^3) ’ (1^) 


and hence (15) reduces to 


= yy / / dTidT2m^mjfi{Ti)fj{T2) 

, J Jo Jo 

2N 

X i:[r(r, - T2%+NA'^lVj) . 


(17) 


1=1 


The parameter a is obtained by substituting (14) and (17) in 
Eq.(lO). We have the most general form now for the work 
fluctuation and its mean value in terms of the dynamics of the 
system. The Eqs.(9), (14) and (17) are the key results of this 
paper. 


III. CONDITIONS FOR a TO BE UNIVERSAL CONSTANT 

The parameter a as determined by Eqs.(14) and (17) can 
depend on the details of the system and the forces. Clearly 
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we need to find conditions when is proportional to {W). 
From Eq.(17) it is clear that we need at least the steady state 
fluctuations such that 

{ipiVj) = Dj5ij+N , Dj = (vj), (18) 

then tr^ reduces to 

/ / dTidT2MTi)fj{T2) 

. ^ Jo Jo 

X Vt+Nj+N{Ti - T2)Djmimj . (19) 

If the force acts only on one variable say, i = k i.e. /j = 5i^f, 
then (14) and (19) reduce to 

cr^ = 2mD^{W) , (20) 

2{W) = f f dTidT2f^{Tl)f^{T2) 

Jo Jo 

X V^+N,K+N{Ti-T2)m. (21) 

In this case a reduces to 

a = 2{W)/2{W)D^ = 1/mD, . (22) 

The equipartition theorem gives mDn = KT and hence 

a=l/KT. (23) 

Another possibility is when rrijDj is independent of the index 
j, then a reduces to the universal constant. Thus the condi¬ 
tions when a reduces to a universal constant for a multidimen¬ 
sional system are (I) Eq.(18) and the force acts only on one 
variable or (II) rrijDj must be independent of the index j. Eor 
a system with one degree of freedom like a harmonic oscilla¬ 
tor, like that obtained by linearizing the non-linear potential 
around some mean value, we have a steady state described by 

exp{— -I- /KT}. Then the steady state mean 

value is such that (xp) = 0. The above conditions are satisfied 
and the ET with ol = 1/KT holds almost trivially. 


and 


= dj - VkSkj - 


(26) 


where j = x, y. Here 7 is the damping coefficient, ujc is the 
cyclotron frequency, 5 = ^ ^ ^. The pj ’s are Gaussian 

processess with 


hliit)) = 0 (27) 

{pjit)pji{t')) = 5jji5{t — t')2m'yKT. (28) 

In the steady state the distribution is proportional to 

exp I - i^rn[vl + v^] + + y^])/KT | (29) 

and hence 

{vD = {vD = = n^y^) = KT/m. (30) 

All terms like {vxVy) etc are zero. Let us now consider a force 
on the electron say by a time dependent electric field in x di¬ 
rection. Thus the work variable is 

w = m [ drf{r)Vj^{r) . (31) 

The Eqs. (14) and (17) for the present problem reduce to 

(iV = 2) 


{W)=mf dri f dT2f{Ti)f{T2)V33{Ti-T2), (32) 

Jo Jo 


and 


2 2 
(j^ = m 


[ [ dridT2/(ri)/(r2) - T2)]3,/(^/’zt^i) 

Jo Jo 

[ f dTldT2f{Tl)f{T2)V33{Tl-T2){vl). (33) 
Jo Jo 


IV. FLUCTUATION THEOREM FOR DISSIPATIVE 
CLASSICAL CYCLOTRON MOTION 


Thus 


cr^ = 2m{vl).{W) = 2KT{W) (34) 


We will now discuss the applicability of the general results 
of Sec.II to the case of an electron which is in a magnetic field 
B applied along the z axis. The electron is further confined in 
a parabolic potential[ 12] . The confining potential is especially 
important in the quantum dynamics. The Hamiltonian is given 
by 

^ i (p-W + lrnn\x^+y% (24) 

with Ax = ^By, Ay = —^Bx. Here H is the frequency 
of the trap potential. The Langevin equations for this system 
have been derived earlier and are given by [13] 

Xj = Vj , Vj = (^pj - Aj'^ jm , (25) 


which on substituting in ( 10 ) yields the universal value of 
a = 1/KT. We have thus shown that the FT in the standard 
form is valid for the dissipative cyclotron motion though it is 
described by a four dimensional Gaussian Markov process. 

The matrix element V 33 (r) is obtained from the solution of 
(26) for Vx{t) under the initial conditions Xi ( 0 ) = X 2 ( 0 ) = 0 , 
= 0’ 'fx(O) = 1, p —?■ 0. The Laplace transform of 
^33 (''■) is found to be 

^33(^) = f z + 7 J -1- - —• (35) 

V z z + j+^J 

For H —>• 0, Vsoiz) yields 

V33{t) = 


cos UJcT 


(36) 
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Given the applied force’s functional form, one can compute 
the explicit results for {W) and {W'^), as shown below. 

For simplicity we present below the explicit results for 
n = 0. Neglecting transient terms (proportional to the ini¬ 
tial velocities), the solution of Eq.(26) for Vx{t) reads: 

Vx{t) = - / dTe^'"<[rix{T)cosujc{t-T) 

m Jo { 

- r]y{T) sin OJc{t - t)] 

+ f{t)cosujc(,t - t) \ . (37) 


Thus, 


{vx{t)) = - - / cos Wc(i-r), (38) 

m Jo 

and therefore, 

(W{t)) = — [ dti f 
m Jo Jo 

X cosujciti — h) ■ (39) 

Our next query is with regard to (t)) which is given by 

dti [ dt2f{ti)f{t2){Vx{ti)Vx{t2)), (40) 
Jo 

which, for the sake of convenient time ordering (G > t 2 ), can 
be rewritten as: 

dti [ dt2f{tl)f{t2){Vx{tl)Vx{t2), 

Jo 

(41) 

where the velocity-velocity correlation is to be obtained from 
Eq.(26). In the steady state (G, t 2 —>■ oo) with G — f 2 = finite, 
we find 

9KT /■* 

{W^t)) = - dh df2/(ii)/(f2)e-^(‘^-*^) 

m Jo Jo 

X cos Wc(fi — ^ 2 ) • (42) 

Comparing with Eq.(39), we find 

=2KT{Wit)), (43) 


(W^{t)) = 2 f 
Jo 



and hence, the parameter a in Eq.(lO) is identical to: 


a = {KT)-^. 


(44) 


Thus, 


p{-W)/p{W) = exp(-VF/fcr), (45) 


validating the GCET, in the present case. 

V. CONCLUSIONS 


We have derived a class of generalized fluctuation theorems 
valid for a system described by a Gaussian Markov process. 
The derivation is fairly broad-based and thus this class of ETs 
would also be valid for non-equilibrium steady states and not 
just for thermal equilibrium. We then present conditions on 
the Markovian dynamics so that we recover standard form of 
the ET. We apply the results of Secs II and III to dissipative 
cyclotron motion and show the applicability of the ET even 
though the dissipative dynamics is in a four dimensional phase 
space. Our general results would be required, in the case of 
the so-called “isometric” fluctuations, if the trapping potential 
is anisotropic and the friction coefficients in x and y directions 
are non-identical[9-l 1]. 

To this date much of the applications of the ET have been to 
classical, non-equilibrium and driven systems in the realm of 
complex fluids and biological processes[14-18]. In as much 
as many of these systems are governed by Brownian and gen¬ 
eralized Brownian dynamics, the present results for general¬ 
ized Gaussian Markov processes are expected to have useful 
relevance. 

In a subsequent paper, we shall consider the extension of 
the present study to a hitherto less-explored area of quantum 
fluctuations, with specific focus on the dissipative cyclotron 
motion. Most of the classical results derived for the latter case 
would emerge as natural, special expressions in the appropri¬ 
ate classical limit. 

We thank Ajay Sood for useful discussions. GSA thanks 
M. Barma, Director, TIER for hospitality during the course of 
this work. 
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